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Abstract — We construct families of high performance quantum 
amplitude damping codes. All of our codes are nonadditive and 
most modestly outperform the best possible additive codes in 
terms of encoded dimension. One family is built from nonlinear 
error-correcting codes for classical asymmetric channels, with 
which we systematically construct quantum amplitude damping 
codes with parameters better than any prior construction known 
for any block length n > 8 except n = 2 r — 1. We generalize this 
construction to employ classical codes over GF(3) with which 
we numerically obtain better performing codes up to length 14. 
Because the resulting codes are of the codeword stabilized (CWS) 
type, easy encoding and decoding circuits are available. 



I. Introduction 

Quantum computers offer the potential to solve certain 
classes of problems that appear to be intractable on a classical 
machine. For example, they allow for efficient prime factor- 
ization [1], breaking modern public -key cryptography systems 
based on the assumption that factorization is hard. Quantum 
computers may also be useful for simulating quantum systems 
[2], [3]. 

However, quantum computers are particularly subject to the 
deleterious effects of noise and decoherence. It was thought, 
for a time, that quantum error-correction would be precluded 
by the no cloning theorem [4] which seems to rule out redun- 
dancy as usually employed in error correction. The discovery 
of quantum error-correcting codes [5], [6] that allow for fault- 
tolerant quantum computing [7] significantly bolstered the 
hopes of building practical quantum computers. 

For the most part, people have concentrated on dealing with 
the worst case — arbitrary (though hopefully small) noise. This 
turns out to be equivalent to correcting Pauli-type errors, a x = 

(l D'Vy = (i ~o)> a z = (o -l)' actin § on a bounded-weight 
subset of the qubits in the code. Since the Pauli operators form 
a basis of 2 x 2 matrices, a code that can correct all Pauli errors 
can in also protect against any general qubit noise [8], [9]. 

However, as first demonstrated by Leung et al. [10], design- 
ing a code for a particular type of noise can result in codes 
with better performance. In practice the types of noise seen 
are likely to be unbalanced between amplitude (er^-type) errors 
and phase (<7 2 -type) errors, and recently a lot of attention has 
been put into designing codes for this situation and in studying 
their fault tolerance properties [11] [12] [13] [14]. 
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In this paper, we will focus on amplitude damping noise, 
another type of noise seen in realistic settings. Amplitude 
damping noise is asymmetric, with some chance of turning 
a spin up |1) qubit into a spin down |0) state but never 
transforming |0) to This models, for example, photon loss 
in an optical fiber: A photon in the fiber may leak out or 
absorbed by atoms in the fiber, but to good approximation 
photons do not spontaneously appear in the fiber. Several 
people have considered this type of noise [10], [14], [15] but 
there is no systematic method for constructing such codes. 
In general it is a difficult problem to design codes for any 
particular noise model. 

In this paper we present a method for finding families 
of codes correcting one amplitude-damping error. We begin 
with an ansatz relating a restricted type of amplitude-damping 
code to classical codes for the binary asymmetric (or Z-) 
channel. The Z-channel is the classical channel that takes 1 to 
with some probability, but never vice versa^] The amplitude 
damping channel is its natural quantum generalization. The 
problem of designing codes for the amplitude damping channel 
is thus reduced to a finding classical codes for the Z-channel, 
subject to a constraint. This lets us carry over many known 
results from classical coding theory. 

We further simplify the problem by using a novel mapping 
between binary and ternary codes. This allows us to find 
quantum amplitude-damping codes by studying ternary codes 
on a greatly reduced search space. 

The rest of the paper is organized as follows. In section 
PI we describe quantum channels and the quantum error- 
correction conditions. In section [Til] we define what it means 
to correct amplitude damping errors and show how they relate 



to classical symmetric codes. In section IV we show how 
a particular class of amplitude-damping codes arises from 
classical codes for the asymmetric channel, and give some 
new codes based on powerful extant results on classical Z- 
channel codes [16]. In section [V] we define a mapping from 
binary to ternary codes (and back) and use this to construct 
new and better amplitude damping codes. Finally, in section 
VI we summarize our results and give a table of the best 
amplitude-damping codes and how they compare to previous 
work. 

II. Preliminaries 

Pure quantum states are represented by vectors in a complex 
vector space. We will be concerned with finite-dimensional 

'Not to be confused with quantum <r z errors, the channel takes its name 
from its diagram resembling the letter 'Z.' See Figure [T] 
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systems. The simplest quantum system (called a qubit) can 
be described by an element of C 2 , and n qubits together are 
described by an elements of C 2 <g> . . . (g> C 2 = (C 2 )®". Such 
pure states are always chosen to be normalized to unity. More 
generally a quantum system can be described by a density 
matrix, a trace one linear operator from (C 2 )®" to (C 2 )®", 
usually denoted p. 

The most general physical transformations allowed by the 
quantum mechanics are completely positive, trace preserving 
linear maps which can be represented by the Kraus decompo- 
sition: 

U{p) = A kpA\ where ^ A\A k = 1. (1) 

k k 

For example the the Kraus operators for the depolarizing 
channel, the natural quantum analogue of the binary symmetric 
channel, are the Pauli matrices. The Kraus operators for the 
amplitude damping channel with damping rate e are 

A quantum error correcting code is subspace of (C 2 )® ra 
which is resilient to some set of errors acting on the individual 
qubits such that all states in that subspace can be recovered. 
For a rf-dimensional codespace spanned by the orthonormal 
set \tpi), i — 1 . . . d and a set of errors £ there is a physical 
operation correcting all elements E^ G £ if the error correction 
conditions [17], [8] are satisfied: 

Vij>„ (ipilE^Evlipj) = CftuSij, (3) 

where depends only on fj, and v. 

III. Correcting amplitude damping 

For small e, we would like to correct the leading order errors 
that occur during amplitude damping. Letting A = a x + i<r y , 
B = I — <j z , we have 

Ai = ^A, A = I- £ -(I-a z ) + O(e 2 ). (4) 

It can be shown if we wish to improve fidelity through an 
amplitude damping channel from 1 — e to 1 — e* it is sufficient 
to satisfy the error-detection conditions for 2t A errors and t 
a z errors. We will say the such a code corrects t amplitude 
damping errors since it improves the fidelity, to leading order, 
just as much as a true i-error-correcting code would for the 
same channel. We will use the notation \\n, K, ij J to mean an 
n-qubit code protecting a Ji'-dimensional space and correcting 
t amplitude damping errors, sometimes referring to this as a 
t-AD code. Our notation descends from the traditional coding- 
theory notation of [n, k, d] to mean an n-bit classical code of 
distance d protecting k bits and [[n, k, d]] to mean an n-qubit 
quantum code of distance d protecting k qubits. Note that our 
AD notation uses K as the full dimensions of the protected 
space, not k, the log of the dimension. This is in preparation 
for the codes we will design which do not protected an integral 
number of qubits. 

Since the amplitude damping channel is not a Pauli channel 
the usual tools for designing quantum codes cannot be directly 



used. One possible approach would be to design CSS [5], 
[6], [18] codes with different a x ,a z distances [19]. For the 
particular case of single-error-correcting AD code, we then 
would like to have CSS code of a x distance 3 (correcting a 
single a x error) and a z distance 2 (detecting a single a z error). 
Gottesman gives a construction of this kind of CSS code in 
Chapter 8.7 of [20]. We summarize his result as follows: 

Theorem 1 If there exists a binary [n, k, 3] classical code C 
and 1 (the all 1 string of length n) is in the dual code of C, 
then there exists an |_|_^)2 fe_1 , lJJ code. 

These codes indeed have better performance than codes 
designed for depolarizing channels. For instance, a [[7, 2 3 , lj J 
exists while only [[7,1,3]] single-error-correcting stabilizer 
codes exist for the depolarizing channel. In general, the 
classical Hamming bound for [n, k, 3] codes gives k < n — 
log(n + 1), which gives a bound for [[n, k]] single-error- 
correcting AD codes constructed by Theorem [T] i.e. 

fe<n-l-log(n+l), (5) 

while the quantum Hamming bound (cf. [20]) gives 

k < n-log(3ra + l) (6) 

for [[n, k, 3]] stabilizer codes for the depolarizing channel. 

However, one expects that these codes cannot be optimal; 
since we only need to correct a x +ia y , correcting both a x and 
a y is excessive and would seem to lead to inefficient codes. 
Fletcher et al. took the first step toward making AD codes 
based on the non-Pauli error model, i.e. codes correcting a x + 
io y error, not both a x and a y errors [14]. Their codes are 
stabilizer codes with parameters [[2n, n — 1]] and correct a 
single amplitude damping error. Later another work [21] took a 
further step toward making AD codes correcting a x +ia y error. 
These works constructed some nonadditive codes correcting a 
single amplitude damping error, and via numerical search for 
short block length found AD codes with better performance 
than codes given by the CSS construction of Theorem [T] 

The construction of [21] consists of codewords \ip u ) of the 
self-complementary format [22], which is 

km = -j= i\u) + m , (7) 

where u is a binary string of length n and u = 1 © u. 

As observed in [22], which focused on nonadditive single- 
error-detecting codes, codes consisting of codewords given by 
Eq. (|7| automatically detect a single a z error, so we have, as 
shown in [21]: 

Theorem 2 A self- complementary code corrects a single am- 
plitude damping error if and only if no confusion arises 
assuming the decay occurs at no more than one qubit. 

We will take the above observation as a starting point for 
making amplitude damping codes, by choosing classical self- 
complimentary codes which correct single errors arising from 
the classical asymmetric channel (or Z-channel). 



3 



IV. Systematic construction from classical 

ASYMMETRIC CODES 

Now we would like to relate the self-complementary con- 
struction to classical error correcting codes for the asymmetric 
channel. Before doing that we first briefly review the classical 
theory of those codes. 

Definition 1 The binary asymmetric channel (denoted by Z 
in Fig. |7J is the channel with {0, 1} as input and output 
alphabets, where the crossover 1 — > occurs with positive 
probability p, whereas the crossover 1 — ► never occurs. 




Fig. 1. The binary asymmetric channel Z and the ternary channel T . 

We will call a classical code that protects against one error in 
the binary asymmetric channel Z a l-code and use the notation 
[n, K, t\ analogous to our notation for the quantum amplitude 
damping code. 

We can then formalize our observation as: 

Theorem 3 If C is a classical [n, K, lj code and^u G C, u G 
C, then Q — {\u) + \u), u G C} is a single-error correcting 
amplitude damping code, [[n, K/2, 1JJ. 

This theorem is almost a direct corollary of Theorem [2] so 
we omit a detailed proof. The main idea is that a classical code 
C that contains both u and u takes care of correcting amplitude 
damping errors while the self-complementary form of 
takes care of detecting the phase errors. And the size of the 
quantum code Q is of course K — \C\/2. This theorem allows 
us to use any classical self-complimentary l-code to construct 
self-complementary amplitude damping codes. The question 
that remains is how to find classical self-complimentary 1- 
codes. 

Varshamov showed almost all linear codes that are able to 
correct t asymmetric errors are also able to correct t sym- 
metric errors [23]. Therefore, to go beyond i-symmetric-error 
correcting codes, we will look to non-linear constructions. 
Note that the quantum codes we construct from these non- 
linear codes are codeword stabilized codes, so these nonlinear 
classical codes will typically result in nonadditive quantum 
codes [24]. 

A. Constantin-Rao Codes 

Constantin-Rao (CR) Codes [16] are the best known non- 
linear 1-codes. These beat the best symmetric single-error- 
correcting codes for all n + 1 2 r — 1. An n-bit CR codes is 



constructed based on an abelian group G of size n + 1. The 
group operation is written as '+' for abelian groups. 

Definition 2 The Constantin-Rao code C g \/g G G is given by 

n 

C g = {{(x 1 ,x 2 ,--,x n )\'^2x l gi =g mod n + l}), (8) 
i=i 

where Xi G {0,1} and Qi,g2, •••,<?« are the non-identity 
elements of G. 

The cardinality of C g is lower bounded by 

, , 2" 



for some g G G. 

Let o(g) be the order of g, then it is known 

\C \ > \Cg\, (10) 

with equality if and only if o(g) is a power of 2. 

For a given nonprime n+l, there may be many abelian 
groups of size n + l. If the group G is a cyclic group of order 
n+l, then the corresponding codes are called Varshamov- 
Tenengol'ts codes [25]. It is known that the largest Constantin- 
Rao code of length n is the code Cq based on the group G = 

© P |„+i ©r=i ^ where " + 1 = ^ P \n+iP np [26]. 

An exact expression for the size of a CR code based on the 
group properties is known, and a basic result is that for any 
group G and any group element g, \C„ \ has size approximately 
(for a review, see [26]). Note is the Hamming bound 
for 1-error correcting codes over the binary symmetric chan- 
nel. Thus, CR codes provide excellent performance compared 
to symmetric codes and, indeed, outperform the best known 
symmetric codes for all block-lengths but n = 2 r — 1. 

B. Amplitude damping codes from Constantin-Rao codes 

To build quantum codes from C g , we need to find CR codes 
which are self-complimentary (and preferably large). We will 
show these exist for all n > 1. 

Fact 1 For even n, the Constantin-Rao code Cq is self- 
complementary. 

This is based on a simple observation that all the nonzero 
group elements add up to zero for any abelian group of even 
size. 

The case of odd lengths n is more complicated. We first 
consider the case where n = 4k + 3. Recall that the largest 
Constantin-Rao code of length n is the code Co based on the 
group G = © p |„ +1 -'=i Z p , where N = n p |„ +1 p'V Then 
further note that for an abelian group Z2 © Z2 © G, where the 
group G is of odd size, all the nonzero group elements add 
up to zero. This leads to the following 

Fact 2 For n = 4k + 3, the Constantin-Rao code Cq of the 
maximal cardinality is self-complementary. 

Since |Cq| > \C g \ > ^rj, AD codes constructed from Fact 
[T] and Fact [3] outperform the CSS AD codes of even length 
and odd length n = 4k + 3 constructed by Theorem [T] 
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Note we also have 
Fact 3 For n = 4k +3, the Varshamov-Tenengol'ts code V .-n 

4 

of the maximal cardinality is self-complementary. 

The case for n = 4k + 1 is more tricky. We cannot 
directly get a self-complementary code of length n from some 
Constantin-Rao codes C g of the same length n. But instead we 
can construct self-complementary AD codes of length n from 
the Varshamov-Tenengol'ts codes V g of length n + 1. 

Fact 4 For n = 4fe + l, f/ze shortened Varshamov-Tenengol'ts 
code V n+2 -r obtained by deleting an odd coordinate r from 

2 

Varshamov-Tenengol'ts code V n+2- r of length n + 1 is se/f- 
corap/eraenfary. 

The codewords of this shortened Constantin-Rao code are 
given by 

En + 2 — r 
ixi — mod n + 2. (11) 



we have 



i mod n + 2 



n - 



2 — r, for any set of a^s 



E 



mod n + 2 



(12) 



where Xj <E {0, 1} and ajj = l(BXi. If the a^s satisfy (Hi then 
so do the XiS. Therefore V' n+2 - r is self-complementary. 

2 

It is known that the size of these shortened Varshamov- 
Tenengol'ts codes are approximately [26]. But we know 
that the size of binary symmetric codes for length n = 4k + 1 
is upper bounded by [27], so the construction of AD 
codes given by Fact [4] also outperforms the CSS AD codes of 
length n = 4k + 1 constructed by Theorem [T] 

Example 1 For n = 8, choose the abelian group of size n + 
1 = 9 be Z 3 © Z 3 . The codewords of the Constantin-Rao code 
Cq are given by a linear code C\ generated by 



{00000011, 00001100, 00110000}; 



(13) 



and four pairs Pi (i=1...4): 



Vi = {10100001,10101101}, 
V 2 = {10000110,10110110}, 

r 3 = {01100100,01100111}, 

V A = {00101010,11101010}; (14) 

and all the complements of Ui=l ^* U^- 

The weight distribution of this code is given by (for defi- 
nition of weight distribution, see [28], [29]) A — l;A\ = 
0;A 2 - 1/4; A 3 = 0;A 4 = 9/2; A 5 = 0; A 6 = 9/4; A 7 = 
0; As — 8. Some of them are non-integers, so this code is 
nonadditive. 

The size of the quantum code is 16, so this is a \_[8, 2 4 , 1JJ 
code. Note the CSS AD code constructed by Theorem [7]/or 
n = 8 gives parameters \_[8, 2 3 ,1JJ. And the best single- 
error-correcting stabilizer code for the depolarizing channel 



is [[8, 3, 3]]. Therefore, this nonadditive AD code encodes one 
more logical qubit than the best known stabilizer code with the 
same length and is capable of correcting a single amplitude 
damping error. 

For short block length (< 16), a comparison of the code 
dimensions given by this Constantin-Rao construction with 
other constructions will be listed in Table U in Sec. [Vi] One 
can see that this Constantin-Rao construction outperforms all 
the other constructions apart from the GF(3) construction 
given in Sec. [v] However, since the GF(3) construction is 
not systematic (those codes given by the GF(3) construction 
in Table [I] are found by numerical search), this Constantin- 
Rao construction is the best known systematic construction 
for single-error-correcting AD codes. 

V. The GF(3) construction and the ternarization 

MAP 

We will begin by defining a channel T which acts on a three 
letter alphabet and find ternary codes on this channel. We will 
then show that such codes are related to binary codes for the 
asymmetric channel and since the binary codes will be self- 
complimentary by construction that they will yield quantum 
amplitude damping codes as well. 

A. The ternarization map 

Definition 3 The ternary channel (denoted T in the fig- 
ure) has {0,1,2} as input and output alphabets, where the 
crossovers ->• 0, -> 1, -» 2, 1 -> 0, 1 ->• 1, 2 -> 0, 
and 2 — > 2 all occur with nonzero probability, but 1 — » 2 and 
2 — » 1 never occur. 

We define a map that takes pairs of binary coordinates into 
a single ternary coordinate. There are four possible values of 
binary pairs, and only three ternary coordinates, so it cannot 
be one-to-one. 

Definition 4 The ternarization map & : F% — » F3 is defined 
by: 

6 : {00, 11} -> 0, 01 -»■ 1, 10 -»■ 2. (15) 

This is not a one to one map. So the inverse map needs 
to be specified carefully, that is, a ternary symbol after the 
inverse map gives two binary codewords 00 and 11. 



Definition 5 The map 6 : F 3 — > F2 is defined by: 
& : -> {00,11}, 1 — ^ 01, 2 -> 10. 



(16) 



For a binary code of length n = 2m, by choosing a pairing 
of coordinates, the map 6 m : F| m — > F™ then takes a given 
binary code of length 2m to a ternary code of length m. 



Example 2 The optimal l-code of length n = 4 and 
dimension 4 has four codewords {0000, 1100, 0011, 1111}. By 
pairing coordinates {1, 2} and {3, 4}, the ternary image under 
6 2 is then {00}. 



On the other hand, S 1 



F5 



Fjj takes a given ternary 



code of length m to a binary code of length 2m. 
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Example 3 By starting from the linear ternary code [4, 2, 3] 3, 
with generators {0111, 1012}, we get the binary image code 



C (8) under 6 4 : 








UU (JUL) (JUL) 


nnnnnm 1 
(JUL) (JUL) 11 


nnnm 1 nn 


1)1)1)1)1111 


00110000 


OOllOOll 


00111100 


OOllllll 


11000000 


11000011 


11001100 


llOOllll 


11110000 


11110011 


11111100 


liiinii 


00010101 


00101010 


11010101 


moioio 


01000110 


10001001 


01110110 


10111001 


01011000 


10100100 


01011011 


10100m 


10010001 


01100010 


10011101 


OllOlllO 



(17) 



which is of dimension 32 and corrects one asymmetric error. 
Note this gives exactly the same binary l-code as the one given 
in Example [7] which is the Constantin-Rao code Cq of length 
n = 8 constructed from the group r L-j,® r L^. This example hints 
at some relationship between the GF(3) construction and the 
Constantin-Rao codes. 

B. The GF(3) construction for l-codes 

1) Even block length: Example [3] suggests that good 1- 
codes may be obtained from ternary codes under the map 6 m . 
We would like to know the general conditions under which a 
ternary code gives a l-code via the map & n . The main result 
of this section states that any single-error-correcting code for 
the ternary channel T gives a l-code under the map 6 m [30]. 

It will be useful in what follows to define an asymmetric 
distance between two codewords: 

Definition 6 Letting N(x,y) = = and yi = 1}, we 

define the asymmetric distance between x and y as 

A(x,y) :=max{iV(x,y),iV(y,x)}. (18) 

It is easy to see that a set of codewords with minimum 
asymmetric distance 2 is a l-code. 

Theorem 4 If C is a single-error-correcting ternary code for 
the channel T of length m, then C — & m (C') is a l-code of 
length 2m. 

Proof For any two ternary codewords c'^c?, G C, we need 
to show that the asymmetric distance between & m (c' l ) and 
S m (c 2 ) is at least two. 

First, we cover the case when c[ = c' 2 . Distinct binary 
codewords may arise from the same ternary codeword due to 
the two different actions of & on 0. Such codewords have 
A > 2 since A(00,11) = 2. 

Next, if the Hamming distance between and c 2 is 
three, then the distance between ©"'(c^) and 6 m (c 2 ) is 
also three since A(00, 01), A(ll, 01), A(00, 10), A(00, 01), 
and A (01, 10) are all one and three such As occur. 

Finally, the following Hamming distance two pairs are 
allowed in a single-error-correcting ternary code for T: 

01,22 10,22 01,12 10,21 02,11 

20,11 02,21 20,12 11,22 12,21 ( ' 

It is straightforward to verify that 6 on these pairs also results 
in binary codes with A > 2. □ 



The following corollary is straightforward. 

Coroliary 1 If C is a linear [n, fc,3]3 code (the subscript 
indicates that the code is over a three-letter alphabet rather 
than a binary alphabet), then & m (C) is a l-code of length 
2m. 

2 ) Odd block length : Theorem [4] only works for designing 
l-codes of even length. Now we generalize this construction 
to the odd length situation, starting from 'adding a bit' to the 
ternary code [30]. 

Definition 7 We call a code acting on F 2 x F™ a generalized 
ternary code of length m+1. We further adopt the conventions 
that & m (C') gives a (2m+l)-bit binary code by acting on the 
m trits of a generalized ternary code C and & 2m (C) when C 
has length 2m + 1 gives a generalized ternary code by acting 
on the last 2m bits of C. 

Theorem 5 If C is a single-error-correcting generalized 
ternary code for the channel Z x T rn of length m + 1, then 
C = & m (C') is a l-code of length 2m + 1. 

Proof 

As in the proof of Theorem |4] we need to show that for 
any two codewords c' l7 c 2 6 C, we need to show that the 
asymmetric distance between ©"(c^) and 6 m (c' 2 ) is at least 
two. If the Hamming distance between codewords on just the 
ternary part of the code is at least two, then the situation 
reduces to the previous proof. 

We need only worry about the case where the Hamming 
distance between c[ and c 2 is two, and one of the differences 
in on the binary coordinate. Assume the first coordinate is a 
bit and the second is a trit, then since C' is a single-error- 
correcting generalized ternary code the only allowed pairs 
are 01,12; and 12,11. The corresponding images of each 
pair under 6™ give binary codewords of asymmetric distance 
A = 2. □ 

To illustrate this generalized ternary construction, let us look 
at the following example. 

Example 4 The code {0000,0111,0222,1012,1120,1201} 
corrects a single error from the channel Z x T 3 . Under the 
map 6 3 it gives the binary code 

0000000 0000011 0001100 0001111 
0110000 0110011 0111100 0111111 

ooioioi oioioio looono mono ' 

1011000 1011011 1100001 1101101 

which is a binary code of length 7, dimension 16 which 
corrects one asymmetric error. 

The following corollary is straightforward, but gives the 
most general situation of the ternary construction. 

Corollary 2 If C is a ternary single error correcting code of 
channel Z mi x T'™ 2 of length mi + m 2 , then C = 6" 12 (C) 
is a l-code of length mi + 2m 2 . 
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C. The Gi 7 '(3) construction for AD codes 

1) Even block length: We first examine under which con- 
ditions the image of a ternary code under & could be self- 
complementary. 

Definition 8 A ternary code C is self-complementary if for 
any c' <E C , c' G C , where c = (30c) mod 3 (3 = 33 ... 3, 
i.e. the all '3' string). 

Example 5 The ternary code C = {000,111,222} is self- 
complementary. For 111 G C, TTT = 333 111 = 222. 

Definition 9 We say that binary code C of even length n = 
2m has ternary form if 6 m (S m (C)) = C. 

The properties of S gives the following 

Fact 5 If a ternary code C of length m is self- complementary, 
then its binary image under &, C = S m (C), is self- 
complementary. On the other hand, if a binary code C of 
length 2m is of ternary form and is self-complementary, then 
its ternary image 6 2m (C) is self-complementary. 

To use Fact [5] to construct good single-error-correcting AD 
codes for even block length, first recall Example [T] (and 
Example [3j: 

Example 6 The code given in Example [7] under the & map 
(pairing up coordinates {1, 2}, {3, 4}, {5, 6}, {7, 8}) gives a 
linear code over GF(3) generated by {0111, 1012}. 

We know that all the linear ternary codes are self- 
complementary, so the 1-codes constructed from linear ternary 
codes of distance 3 can directly used to construct single-error- 
correcting AD codes [30]. Since in general we search for self- 
complementary ternary codes C with largest possible size of 
C = S(C'), those AD codes obtained from linear ternary codes 
of distance 3 are sub-optimal. 

We now show that the AD codes given by the Constantin- 
Rao construction are actually a special case of the GF{3) 
construction. 

Theorem 6 For n even, the Varshamov-Tenengol'ts code Vo, 
and the Constantin-Rao code Co of largest cardinality has 
ternary form. 

Proof We only need to prove that there exists a choice of 
pairing, such that for any codeword v G Vo (Co), if v restricts 
on one chosen pair a is 00, then there exists another codeword 
v' G Vo (Co) such that v' — v\g, and v'\ a — 11. Here a denotes 
all the other coordinates apart from a. 

For the Varshamov-Tenengol'ts code Vo of even length n, 
choose the pairing {i, n — i + l}"f J, then the above condition 
is satisfied. This is because i+n — i+1 = n+1 mod n+1 = 0. 

For the Constantin-Rao code Co of largest cardinality, 
which is given by the group G = © r ©7=1 Z Pr ., note n 
is even, so n+1 is odd. Therefore all p r are odd for 
p r \n + 1, where n + 1 = H Pr \ n+ ip" r . Write any group 
element as (sn, Si ni , S21, S2n 2 — ). Then we can pair 



it with {pi - sn, ...,|>i - s lni ,p 2 - s 2 i,...,p 2 - s 2 „ 2 ...), 
mod {pi,...,pi,p2,...,P2,—), where s rjr G {0, ...,p r - 1} 
and j r = 1, n r . □ 

From both Fact Q] and Theorem [6] we learn that for even 
block length, the Constantin-Rao code Co of maximal cardi- 
nality is both self-complementary and has ternary form. There- 
fore, the AD codes given by the Constantin-Rao construction 
is actually a special case of the GF{3) construction. 

2 ) Odd block length: For n odd, we need to generalize the 



GF(3) construction. As already discussed in Sec. V-B.2 for 
n = 2m + 1, we design codes correcting a single error of 
the channel Z x T m . And we call these codes 'generalized 
ternary.' 

We need to examine under which condition the image of a 
generalized ternary code under & is self-complementary. 

Definition 10 A generalized ternary code C of length 2m + 1 
is self-complementary if for any c' G C, c' G C. Here c! x 
I © c[, c' t = 3 e c- mod 3, for i = 2, . . . , m + 1. 

Example 7 The generalized ternary code C = 
{000,100,011,122} is self-complementary, because 
000 = 100 and 0TT = 122. 

The properties of & give the following: 

Fact 6 If a generalized ternary code C of length m + 1 
is self- complementary, then its binary image under the map 
C = S m (C) is self-complementary. On the other hand, if 
a binary code C of length 2m + 1 has generalized ternary 
form and is self-complementary, then its image 6 2m (C) is 
self-complementary. 

We now show that the AD codes given by the Constantin- 
Rao construction are actually a special case of the generalized 
ternary construction. 

Definition 11 A binary code C of odd length n = 2m + 1 has 
generalized ternary form if & m (& m (C)) = C. 

Based on this definition, if a binary code C of odd length 
2m+ 1 has generalized ternary form, then it can be constructed 
from some codes correcting a single error of the channel 
Z x T m via the ternarization map. The following theorem 
then shows that certain Varshamov-Tenengol'ts-Constantin- 
Rao codes are a special case of asymmetric codes constructed 
from single-error-correcting codes for the channel Z x T m 
[30]. 

Theorem 7 For n odd, the Varshamov-Tenengol'ts code V g 
has generalized ternary form. 

Proof We only need to prove that there exists a choice of 
pairing which leaves a single coordinate as a bit, such that for 
any codeword v G V g , if v contains the paired bits 00, then 
there exist another codeword v' G V g which is identical except 
that the 00 pair is replaced by 11, and vice versa. 

For the Varshamov-Tenengol'ts code V g of odd length 
choose the pairing {i, n — i 



IjlHi 1 ^ 2 ' leave the coordinate 
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(n + l)/2 as a bit, then the above pairing condition is satisfied. 
This is because i + (n — i) + 1 = (n+ 1) mod [n + 1) =0. 
□ 

Now recall Fact [5] which states that for block length n — 
Ak + 3, Vn+i is self-complementary. We further show the 
following: 

Fact 7 For n — Ak + 3, Vn+i is of generalized ternary form. 

4 

To see this, do the pairing {i, n — i + ^ . Here we leave 
the coordinate (n + l)/2 unpaired so it is unchanged under 
the map & m . 

For length Ak + 1, recall Fact |4] that the shortened 
Varshamov-Tenengol'ts code V' n+2 - r obtained by deleting any 

2 

'odd' coordinate r from Varshamov-Tenengol'ts code V n+2-r 

2 

of length n + 1 is self-complementary. We further show the 
following: 

Fact 8 For n = Ak + I, the shortened Varshamov-Tenengol'ts 
code V' n+2 -r obtained by deleting any 'odd' coordinate r 
from Varshamov-Tenengol'ts code V n+2- r of length n+l has 
generalized ternary form. 

To see this, for the shortened Varshamov-Tenengol'ts code 
given by 

n+2 n + 2- r 

ixi = mod n + 2, (21) 

i— l,i^r 

do the pairing {i, n — i + 2}™f \. Here we leave the coordinate 
n — r + 2 unpaired so it is unchanged under the map 6 m . 

VI. Summary of new constructions for amplitude 

DAMPING CODES 

For short block length we summarize the results of single- 
error-correcting AD codes obtained from the GF(3) construc- 
tion in Table [I] and compare them with AD codes obtained 
from other constructions. 

Note the [[12,168, lJJ code in Table [i] is cyclic, which 
can be obtained by the classical 1-code [12, 336, lj given in 
[30]. The L Ll , 49 , 1 J J code is 'almost cyclic', from which 
(deleting 4 classical codewords then add another 2) we can 
obtain a cyclic code [[10, 47, lJJ, with classical codewords 

00000 11111 22222 21100 20111 (22) 

and their cyclic shift, plus all the complements. There is 
another cyclic code ((10,47)), with classical codewords 

00000 11111 22222 21100 21011 (23) 

and their cyclic shift, plus all the complements. 

Table [I] shows that the Constantin-Rao construction C g 
outperforms other constructions apart from the (generalized) 
GF(3) construction. This is reasonable since we know that 
the Constantin-Rao construction is actually a special case of 
the (generalized) GF(3) construction. For all lengths up to 
14, the (generalized) GF(3) construction indeed gives AD 
codes of best parameters. Lengths > 14 are out of reach of 



TABLE I 

Codes: This table compares the various constructions for 
amplitude damping codes, giving the best known codes created 
by various constructions. the first column gives the number 
of qubits. the second column gives additive codes. the third 
column uses the construction given in gottesman [20]. the 

third column gives codes created by the complementary 
construction of lang and shor [21]. the fourth column (c 3 ) 
gives constantin-rao codes. the fifth column gives codes 
constructed using theoremfjjand computer search. 



n 


GF(4) 


[20] 


[21] 


c g 


GF(3) 


4 


1 


1 


2 


2 


2 


5 


2 


2 


2 


2 


2 


6 


2 


4 


5 


5 


5 


7 


2 


8 


8 


8 


8 


8 


8 


8 


12 


16 


16 


9 


8 


16 


18 


23 


21 


10 


16 


32 


11 


17 


19 


11 


32 


64 


78 


86 


89 


12 


64 


128 


146 


158 


168 


13 


128 


256 


273 


274 


291 


14 


256 


512 


515 


548 


572 


15 


512 


1024 


931 


1024 


* 


16 


1024 


1024 


1716 


1928 


* 



the current computational power we have. As we know that the 
Constantin-Rao construction outperform the CSS construction 
for all lengths except n = 2 r — 1, where the binary Hamming 
codes are 'good', it is very much desired to know whether 
the (generalized) GF(3) construction can give us something 
outperforms the CSS construction for the length n — 2 r — 1. 
From [30] we know this is possible for classical 1-codes, but 
it remains a mystery for the quantum case, which we leave 
for future investigation. 

Finally, numerical search also found a [[9,26, lJJ single- 
error-correcting AD code (exhaustively found to be optimal 
among all the self-complementary codes), which cannot be 
obtained from any of the above constructions. Also we have 
found, via random search, a [[10,51,1JJ code, which also 
cannot be obtained from any of the above constructions. 
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